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Abstract 

OhI Just like AdS spacetimes, Lifshitz spacetimes require counterterms in order to 

(— I ' make the on-shell value of the bulk action finite. We study these counterterms using 
the Hamilton-Jacobi method. Rather than imposing boundary conditions from the 

p^ \ start, we will derive suitable boundary conditions by requiring that divergences can 

^ ' be canceled using only local counterterms. We will demonstrate in examples that this 

Jy^ ■ procedure indeed leads to a finite bulk action while at the same time it determines 

ly-v \ the asymptotic behavior of the fields. This puts more substance to the belief that 

lO ' Lifshitz spacetimes are dual to well-behaved field theories. As a byproduct, we will 

t~^ ' find the analogue of the conformal anomaly for Lifshitz spacetimes. 
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1 Introduction 



Lifshitz spacetimes were originally introduced as possible holographic dual descriptions of 
non-relativistic field theories [1, 2] and have since appeared in many different setups, for 
example as IR geometries [3]. Moreover, they have appeared as solutions of string theory 
[4-10] and although they are not yet at the same footing as ordinary AdS spacetimes, it 
is worthwhile to explore to what extend the usual AdS/CFT techniques can be applied to 
Lifshitz spacetimes as well. 

Certain features of Lifshitz spacetimes that have been and still are confusing are its global 
causal structure, the absence of a version of "global Lifshitz" , the nature of the bound- 
ary conditions on the metric and other fields, and indications coming for example from 
Schrodinger holography [11] that one needs non-local counterterms to remove divergences 
in the on-shell value of the action. 

Motivated by this we decided to explore the nature of the divergences that appear in 
Lifshitz spacetimes when computing the on-shell value of the effective action using the 
Hamilton- Jacobi method, which turns out to be more efficient in this case than using the 
Fefferman- Graham expansion, which rapidly becomes quite intractable. 

Normally, in order to perform holographic renormalization, one needs to first say something 
about the boundary conditions for the fields. We will, however, follow a different and novel 
approach. As we will show, if we require that all divergences should be canceled by local 
counterterms, this will automatically enforce particular boundary conditions for the fields. 
More precisely, we will find that particular local covariant quantities made out of the bulk 
fields have to scale in a specific way as we approach the boundary of Lifshitz. With this 
approach, we will also show that for a class of bulk Lagrangians all power law divergences 
can indeed be canceled using only local counterterms. This strongly suggests that Lifshitz 
spacetimes are dual to field theories with a well-defined UV completion. 

In addition, certain ambiguities that appear in the analysis of the counterterms have a 
natural interpretation in the dual field theory in terms of marginal deformations, exactly 
as was the case for AdS/CFT. 

Along the way, we will show that counterterms that had been previously proposed in [12] 
are insufficient to cancel divergences beyond the leading order, and find the analogue of 
the conformal anomaly for Lifshitz spacetimes. 

The outline of this paper is as follows. In section 2 we review the Hamilton- Jacobi method 
and apply it to the non-derivative terms in the boundary effective action. All power-law 
divergent terms in the effective action can be canceled using local counterterms. Sometimes, 
logarithmically divergent terms appear which cannot be canceled using local counterterms, 
and it is precisely these that are responsible for the analogue of the conformal anomaly. We 



also describe the relation between ambiguities that appear and the existence of marginal 
deformations. 

In section 3 we perform a non-trivial consistency check by explicitly computing the on-shell 
action for scalar perturbations of the metric and gauge field to second order. We will find 
that with our counterterms the on-shell action is indeed rendered finite. 

Various subtleties, such as the presence of logarithmic divergences, qualitative dependence 
of the answers on the value of the so-called dynamical exponent z, and issues related to 
the boundary conditions are discussed in the conclusions. 

The appendices contain some background material and a brief description of the extension 
of our methods to the terms containing derivatives. 

Note added: As we were preparing this paper for submission to the ArXiv, the paper 
[13] appeared, which reaches similar conclusions as we do though using different methods. 



2 Holographic Renormalization 



In this section we set up the general framework for computing the counterterm action. We 
begin with a brief review of Lifshitz spacetime and the specific bulk action we shall use. 
After that, we describe the Hamilton-Jacobi method of holographic renormalization and we 
introduce the 'Lifshitz scaling anomaly'. Finally, we explicitly compute the counterterms 
at the level of no spacetime derivatives and contributions to the Lifshitz scaling anomaly. 



2.1 Lifshitz spacetime and the Einstein— Proca action 

Lifshitz spacetime is a proposed gravitational dual to a field theory at a UV fixed point 
with anisotropic (Lifshitz-like) scaling symmetry. 

The configuration of [d^ l)-dimensional Lifshitz spacetime [1] that we consider consists of 
the following metric and vector [14], 



ds' = dr'-e''''dt' + e''dx\ A = y/^cToe'^dt. (2.2) 

This metric is invariant under the so-called Lifshitz algebra [15], which consists of time 
translations, spatial translations, spatial rotations, and anisotropic scaling invariance (2.1) 
(with a simultaneous shift in the radial coordinate r \-^ r — log A). Unlike so-called 



Schrodinger spacetimes, the Lifshitz spacetime is not invariant under Galilean boosts 
X I— ;■ x+vt. We will eventually work in 3+1 bulk spacetime dimensions {d = 3), but we keep 
d arbitrary for as long as possible. These fields comprise a solution to the Einstein-Proca 
action S = Sgrav + Sa, with 

^grav = f d'^+^xy^ {R - 2A) + f d'^^^f^ 2K, (2.3) 

Sa = J d'^+'x^^ (-iFt.^F'"' - Irn^A^A^ , (2-4) 

where we used the convention IGvrG = 1. It should be noted that we could also have chosen 
a different action that has the Lifshitz metric as a solution, see e.g. [1, 14]. We have chosen 
the Einstein-Maxwell theory for its relative simplicity. In order to find Lifshitz spacetime 
as a solution, we must pick our parameters to be 

A = --{z^ + z + A), m^ = 2z, a^ = -2^^. (2.5) 

The equations of motion are 

R^v - \r9^v + t^g^v = \t^u. V^F^^ = mM^ (2.6) 

where T^" = ^H=— ^ is the Proca stress tensor. We wish to add a scalar at some point, 
so let us give the scalar action as well, 
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S^ = J d^'^'xV^ [-^^d.cl^d^cl^ - V^(0)J , (2.7) 

with a potential V^(0) = |/i^ 0^ + f30'^ + f 40^ + ..., which is presumed to be known. It is in 
principle possible to consider different setups, e.g. with direct couplings between the scalar 
field and the vector field such as ip'^A^, but for simplicity we will restrict our attention to 
the simple case (2.7). 

Our aim is to construct a finite on-shell action for an appropriate class of asymptotically 
Lifshitz spacetimes, where this notion will be made more precise in the following. As 
usual, the action given in (2.3) and (2.4) diverges on-shell, and it is necessary to introduce 
a set of counterterms to remove these divergences. In this paper we will assume that 
such counterterms are local in the fields, and we will use the Hamilton- Jacobi method to 
determine their form. We expect the on-shell action to be of the form 

5*01 = 5*100 + r, (2.8) 

where Sioc contains the local power-law divergent terms and F diverges at most logarith- 
mically. We will determine 5*100 by imposing the Hamilton constraint H = 0, where H 



is the radial Hamiltonian corresponding to the Einstein-Proca action, which is derived in 
appendix A; it is given by 

H= f d'^xy/^ {Nn + N^-Ha) , (2.9) 

where S^ is a hypersurface of large but constant r, while N and A^" are the usual lapse 
and shift functions. The momentum constraint is given by 

-Ha = -2D\ab - AaDbE' + FabE^ + nda4> = 0, (2.10) 

where the quantities n"'^, E'^ and tt^ are the canonical momenta dual to the induced metric 
7a6, induced vector Aa and the scalar respectively. The Hamiltonian constraint is 

^ = - [^ab^""' - ^^') - Ie^E. - \l - ^,{DaEr -^ = 0, (2.11) 

where C = R — 2A — jFabF"-^ — \m?AaA"- — \da(j)d"'(j) — V{(j)) is the Lagrangian restricted 

to Tjr- 



2.2 Hamilton— J acobi equation and the Lifshitz-scaling anomaly 

The Hamilton-Jacobi (HJ) equations of motion for a point particle are H = —dtSd and 
p(t) = dSci/dq, where the on-shell action Sd is the action evaluated on the classical path 
with given initial and final conditions. The first HJ equation simply becomes H = 0, while 
the second one is generalized to 

vr«^(r) = ^l^(r), E^ir) = ^^(r), n,{r) = ^^^(r). (2.12) 

The two HJ equations may be combined into what is known as the Hamilton-Jacobi equa- 
tion, 

i^(7a.,A„0;f^,Jf^,f^)=O. (2.13) 

The HJ equation is a functional PDE for the on-shell action. In principle, this equation 
determines the form of the on-shell action 5*01, but it is far too difficult to solve. Since we are 
interested only in the local part 6*100 we can recast the problem in a more tractable form. 
It is useful to introduce the following notation for the 'kinetic' part of the Hamiltonian 
constraint (2.11): 

{^^f {F, G]^- Uclbd - -^lablc^i ^^ (2.14) 

V d-l J d-fab hcd 

_15F5G_1 6F 6G ^ j^ ^E 5G 



such that the Hamiltonian constraint is simply 

= n = {S,uSci}-C. (2.15) 

The bracket {F, G} is symmetric and bihnear in F and G. Therefore we can use the 
sphtting ^ci = ^loc + r to write: 

= {Sloe, Sloe} -C + 2{Sioe, T} + {T , T} . (2.16) 

We define the "local part" of this expression as "Hioc = {S\oc, S\oc} — C The divergent 
part of "Hioc should vanish by itself because the non-local part shouldn't contain power-law 
divergences. Solving Hiocdiv = determines the divergent terms in Sioc, which will be our 
counterterms. 

We will see that, in the presence of marginal deformations, this procedure possibly leaves 
a finite remainder in "Hioc, which we denote by T-Lrcm and is generically determined unam- 
biguously. The presence of such a remainder is refiected in the r dependence of S\oc via 
the Hamilton- Jacobi equation for the local part: 

drSloc =- [d'^X^/^'Hrem, (2.17) 

where the partial derivative on the left hand side is taken keeping the boundary fields fixed. 
Since the HJ equations tell us that drSd = 0, or drSioc = —dj-T, we immediately derive the 
Lifshitz analogue of the Weyl anomaly [16-19]: 

drT = I (fx^/^Urer^. (2.18) 

On the other hand, using the invariance of F under bulk diffeomorphisms [19], we can 
rewrite the left-hand side of this equation as follows (for simplicity we assume ^u = 0): 

drT = [d'x(iab^ + A,^ + <P^^] (2.19) 



d'^x I 2z^uTr- + ^%Tr- + ^^tTT- + K^—]+--- (2-20) 



ST ^ ST ;j ST , ^sr 

Ht S-fij SAt ^ sip. 

d'x./^ {z{Tt') + (T/) + zMOa) + A^ 0(C^)) + . . . , (2.21) 

where we used the hatted notation for the asymptotic values of the fields, e.g. 7tj = 
lim^^oo e~^'^''7tt, A^ is the leading radial scaling of the scalar field and the dots represent 
subleading contributions. We also used 

(r-) = ^^. m^^^, n) = ^g, (2,22) 

V-7 07a6 V-7 SAt V-7 Scj) 

We thus see that the 'Lifshitz scaling anomaly', denoted Az-, is given by 

A = z{Tt') + {T,') + zMOa) +A^0(O^)= hm e(^+2>-Hrem, (2.23) 

SO there is a possibility that the Lifshitz scaling symmetry is broken. We shall see that 
this happens for some specific values of the scalar mass-squared jj?. 



2.3 Initial conditions 

As we mentioned earlier, one typically chooses an Ansatz that is covariant [17-19], such that 
the momentum constraint is automatically satisfied. Solving the HJ equation thus reduces 
to solving the Hamiltonian constraint. In the present case, the most general covariant 
Ansatz one can take is 



>S'ioc = / d Xy/—^ U{a, (j)) + (derivative terms). (2.24) 

The quantity a = AaA"" is the only scalar that one can construct from the metric and the 
vector field containing no derivatives. At present, our main focus will be on this level of 
no spacetime derivatives; see Appendix B for a discussion of the higher- derivative levels. 

As an aside, notice that the induced metric may be viewed as the metric on d-dimensional 
fiat space with an r-dependent speed of light c^ = e'^^"^^'' [12]. More specifically, 

-fab{r) dx^'dx^ = e^'^{-cldf + dx^). (2.25) 

Therefore, heuristically speaking, it seems quite natural to impose covariance on S^; it is 
only in the limit r — )■ oo that the speed of light becomes infinite. 

Before we move on to solving the HJ equation, we need to establish the leading-order 
behavior of U{a, 0) in an expansion about the Lifshitz background. In fact the HJ equations 
are functional differential equations for Sd, and in order to solve them it is necessary to 
provide the initial conditions, that is the value of Sd for a solution of our choice. Henceforth, 
we set d = 3. Recall that the background value of a is given by Oq = —2{z — l)/z, cf.(2.5). 
We shall use the Hamilton equation of the type q = dH/dp in order to fix the leading-order 
behavior of U{a, </>). 

dr'Jab = ^, , „.r = -27rab + 7afe7r, (2.26) 



drAa = ,, ;::i^^, = -Ea + \DaDbE\ (2.27) 

dr(p = ., Zl , = -vr^. (2.28) 

At the level of no spacetime derivatives, the canonical momenta are given (via the second 
HJ equation (2.12)) by: 

1 Ox T 

T^ab = --fabU-AaAb—, (2.29) 

dU 
E, = 2A,— , (2.30) 

dU 
-* = 1^. (2-31) 



K^ 


/-7 7r"'') 
6H 


5(> 


6H 



where now vr'^* = —^-^-^^ denotes only the local part of the momentum (and similarly for 
E"- and tt^), so that we find 

dr7ab = (^U-a ^^ 7a6 + 2A,A ^, (2.32) 

dU 

drAa = -2Aa—, (2.33) 

da 

dU 
M = -^. (2.34) 

Let us expand the function U{a, 0) around the Lifshitz background as 

t/(a,0) = ^M^,(a-ao)'"0". (2.35) 



m,n 



We can then translate the radial behavior of the Lifshitz background dr'~fu = '^^Itt, drlij = 
'^liji drAf = zAt, and dr(f) = (0 = in the background) into the values for the first 
few coefficients in the function U{a,(j)). We evaluate the above Hamilton equation on the 
Lifshitz background, so that we find 

uoo = 2{z + 1), uio = --, uoi = 0. (2.36) 

Therefore, the first three coefficients in the expansion for U{a, (p) around the background 
have been determined by imposing that the background is a solution to the equations of 
motion. 



2.4 Radial behavior of (a — ao) and 

The radial behavior of [a — ao) and </> is dictated by the value of the coefficients U20 and 
U02 respectively. We will now show how this comes about. 

Let us start with the more familiar example of the scalar field 0. The radial behavior can 
be obtained directly from the free scalar field equation in the Lifshitz background, 

d^,(j) +{z + 2)dr<p + 7"' dadbcj) = /iV- (2.37) 

Taking the Ansatz = x(t, x, y) exp(A0 r) yields 

\^ = -^((^ + 2)±v/(^ + 2)2 + v). (2.38) 

Alternatively, one can find this radial behavior by expanding the Hamilton equation (2.34) 
to second order, such that 

dr(f) = -2mo2 + ... (2.39) 



Later, we will find that M02 = — A(^/2, which agrees with the answer provided by the equa- 
tions of motion (2.37). Notice that there is also the Lifshitz analogue of the Breitenlohner- 
Freedman bound, i.e. //^ > —{z + 2)^/4, see also [1, 14]. 

For (a — ao), we similarly find via (2.32) and (2.33) 

dr{a-ao) = --af/ - (a^ + 4a)-— = Xc,{a - ao) + ..., (2.40) 

^ CJCx 

with 

A„ = ^^^^^20-^(^-1). (2.41) 

In this case finding the radial behavior using the Einstein-Proca field equations would be 
considerably more difficult, since a is a composite field containing both the vector and the 
metric. 

One may wonder whether the bilinear term in U{a, 0) spoils these relations. Later on, we 
will find that the HJ equation gives uu = 0, so this term is in fact absent. 



2.5 Non-derivative counterterms 

We shall now set out to find the counterterms at the level of no spacetime derivatives by 
solving the local part of the HJ equation. As we mentioned before, this comes down to 
solving the local part of the Hamiltonian constraint, 

0= {Sloe, Sloe}-/: (2.42) 

-^-^-(r^^-)(f)^4<-Ki)^-'=^-^^'-"-^^'^'' 

(2.43) 

where the symbol = indicates that there might be a finite remainder "Hrem as discussed 
in section 2.2. Remember that S\oc = J df'XyJ—'~^ U + (derivatives) and we expanded the 
function U about the Lifshitz background as f/(a, (p) = ^„^ ^ Umn (a — ao)"*'/'"- We can do 
a similar expansion of the non-derivative part of the local Hamiltonian constraint, 

Hloe,non-deriv = J^ ^-" (" " "o)"*^"" (2-44) 

This allows us to solve the Hamiltonian constraint order by order. 

Order zero and order one. The above values for uqo, mio, and mqi are such that both 
the order-zero and the order-one constraints vanish identically, i.e. "Hoo = "Hoi = "Hio = 0, 
even though they depend on U20 in principle. 
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Order two. The constraints at second order in {a — ao) and/or are 

1 z^ — 1 bz^ 

%o = -2<-(2^-5Ko + 8^^M^o-^, (2.45) 

^11 = \ (s^^Mao - 2mo2 - ^ + 7 j Mil, (2.46) 

z^ -I 1 

^02 = 2 —u\^ + {z + 2)uQ2-2ul^^-^i\ (2.47) 

which is a system of three coupled equations in terms of the three unknowns M20; wh, and 
uo2- There is also one free parameter (apart from z), namely the scalar mass-squared /i^. 
The generic case, for which ^'^ ^ 2{z — l){z — 2) , one finds 



z^ 



U20 



^g^^, _^ (5 - 2. ± V9.2 - 20z + 20) , (2.48) 

nil = 0, (2.49) 

M02 = ^ (^ + 2 ± ^/{z + 2)2 + 4/i2) . (2.50) 

For /i^ = 2(z — 1)(2; — 2), however, one of the three coefficients (^20, ^ii, ^02) remains 
undetermined. In principle, it is possible that there is a remainder in the Hamiltonian 
constraint that can cannot be set to zero by tuning a coefficient. Such a remainder would 
contribute to a Lifshitz scaling anomaly. Although, in the present case, for this specific 
value of /i^ we find no such remainder. 

Remember from Section 2.4 that the coefficients U20 and U02 determine the radial behavior 
for the fields (a — ao) and respectively. The choice of the sign in U20 and M02 correspond 
to choosing the either the normalizable or non-normalizable mode for the fields. For both 
U20 and Uq2, the minus-sign root will correspond to the non- normalized modes. Though, 
for the scalar field there is a special window for /i^ where both modes are normalizable 
and one has the freedom to choose the plus sign in uo2- This special window is given by 
— ^{z + 2)'^ < ^^ < — 1(2 — 2;)(2 + 3z), which is obtained by requiring that the Klein-Gordon 
norm is finite. 
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Order three. The constraints at third order are given by 

^30 = -^(^ + 2-3/3,)m3o, (2.51) 

+ TTT^TZ^ 7^ (35^^ - 160^2 + 2232 - 170 + il^z" - 30^ + 39) /?,) 

128(2"^ — 1)^ ^ ^ ' ' 

H21 = -^(2 + 2-7,-2^,)u2i, (2.52) 

1/ ^ z{z-Q + |3^iz + 2-'y^ , ^ 

-Hu = --(. + 2-/3.-27.)ni2 + ^ ^^TT) ^' ^^'^^^ 

-Hos = -^{z + 2 - 3-f,)uo3 + vs, (2.54) 

where we introduced the abbreviations for the square-roots (3^ = a/Q^^ — 20z + 20 and 
72(/i) = y(F+2p + 4yU^. ^ Recall that v^ is the third order coefficient in the potential 
V{(f)). The ffist thing one sees is that the equations decouple already at third order; this 
continues to hold at higher order (for as long as one does not hit a continuous ambiguity). 
One determines the coefficients M30, M21, M12, and mqs from these equations for generic values 
of n"^. There are two special values for which one of the coefficients remains undetermined. 

• Ambiguity 1. For /i^ = — 1(2; + 2)^, which is above the BF bound, the coefficient 
uo3 remains undetermined. There is a remainder in the Hamiltonian constraint that 
cannot be set to zero, "Hrcm = ^3 0^, there is a contribution to the Lifshitz-scaling 
anomaly (2.23), 

A, = vs lim (e(^+2>^ 03) + . . . (2.55) 

which is present for all values of z. This is the same contribution to the anomaly 
that was found in AdS (z = 1), see e.g. [19, 20] 

• Ambiguity 2. We find a second ambiguity when yU^ = yq{z + 2 — /3^)^ ~" 4(2^ + 2)^, 
the coefficient Uu is undetermined and there is also a remainder, namely 

^'' - M(^Ti) • ^^-^^^ 

Thus, for this specific value of /i^, so we have a contribution to the Lifshitz anomaly, 
A = _±^l±Mi^_tl±M. ,i„ (,<.«..(„ _ „„)^.) + . . . (2.57) 

bA[Z + 1) r^oo 

Interestingly, this contribution to the Lifshitz anomaly vanishes for z = 2. 

^Thc quantity f3z is positive for any value of z. Since 7^ depends on //^, we see that 7^ > as long as 

/.2 >_!(, + 2)2. 
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A„ = ^ ^ > U2o--{z-l), X^ = -2uo2. (2.58) 



a:!: 



2.6 Origin of the continuous ambiguities 

The continuous ambiguities arise when the radial behavior of a term mixing {a — ao) with 
(j) becomes of 0(1). In order to see this, let us recall the radial behavior (at large r) from 
Section 2.4, i.e. {a — ao) ~ exp(Aar) and cj) ~ exp(A<^r), where 

](z^ - 1) 3 

When we plug in the values we found from the HJ equation, we get 

-^ (^ + 2 ± v/(^ + 2)2 + 8(^ -l){z- 2)) , (2.59) 

-^(^ + 2 ± v/(^ + 2)2 + v). (2.60) 

Let us refer to the corresponding modes as (a — ao)± and (f)±. For generic values of 
the scalar mass /i, we see that A~^ describes the radial behavior of the non-normalizable 
modes (sources) and A^ , that of the normalizable modes (vevs). In terms of the scaling 
coefficients, we always have 

X^ + \t = -iz + 2), A^ + A+ = -(^ + 2). (2.61) 

Note that this — (z + 2) cancels against the +{z + 2) coming from a/— 7 ~ e^^^"^^^. 

Order-two ambiguity. The critical value n^ = 2{z — l){z — 2) that we found at second 
order in (a — ao) and (p is simply a result of tuning the scalar mass such that AT = A~. 
We should stress that this ambiguity is not related to the presence of a marginal operator, 
therefore it is not surprising that there is no remainder in this case. This ambiguity is 
parametrized by Uu, which mixes a — Oq and by giving linear contributions to the right 
hand side of (2.39) and (2.40). This suggests the possible presence of a one-parameter 
family of allowed boundary conditions, but we leave this for further study. 

Order-three ambiguities. At third order, we found two ambiguities, which comes from 
the freedom to tune the scalar mass such that either the combination (p^ or (a — ao)0^ is 
of order one in e~^^~^^^'^. Another way of saying this is that either the coefficient Mqs or U12 
remains undetermined. In the first case, i.e. for /x^ = —^iz + 2)^, we indeed see that 

3X- = -{z + 2). (2.62) 

This continuous ambiguity comes from the appearance of a marginal deformation. Namely, 
the operator 0^ is marginal for this specific value of /i^. In the second case, i.e. for //^ = 

^(^ + 2 - (3,y -liz + 2)2, we find that 

A; + 2A^ = -(z + 2), (1<^<2), (2.63) 

A^+2A+ = -{z + 2), {z>2). (2.64) 

In this case, the operator {a — ao)0^ is niarginal. 

13 



3 Renormalized On-Shell Action 



Up to this point we did not impose any boundary conditions on our space of solutions. 
In this section we will analyze the Einstein-Proca equations in the constant perturbation 
sector up to second order and show that divergences in the third-order on-shell action are 
indeed removed using our formalism for a large class of boundary conditions. 



3.1 Perturbative analysis of the Einstein— Proca equations 

First, we discuss the first-order solution obtained in [12] and then we set out to solve the 
second-order equations. The purpose of finding these solutions is to perform a non-trivial 
check of the counterterms that we found in 2.5. Again, we focus on the non- derivative sector 
throughout this paper, so we shall restrict our analysis of the linearized field equations to 
constant modes which only depend on the radial coordinate r. For simplicity, we postpone 
the treatment of vector and tensor perturbations to section 3.3, and here we focus only on 
the scalar subsector. Furthermore, we set the scalar field to its background value as well, 

= 0. 

We adopt the same parametrization as [12] for which the Lifshitz geometry is perturbed 
as follows. 



lu = -e'^'\l + ef{r) + e'fXr) + ...j, (3.1) 

7., = e% (l + 5 A;(r)+52^(r) + ...), (3.2) 

At = V^^oe^'(l + e{j{r) + y{r))+e^j{r) + yXr)) + ...). (3.3) 



We use the small parameter e to keep track of the order in the perturbative expansion. We 
work in radial gauge, which means that the components Qrfj. do not receive any corrections. 

First-order solution. As was noted in [12], the first order field equations for constant 
perturbations reduce to the following three equations. 

= 2f - (2 + 1)/ - {Az + 6)f + 2{z + A){z- l)j, (3.4) 

= (^ + 1)/" + 3{z + l)f' -{z- l)(4z + 2)/ -{z- l){Az^ + 6z + 8)j, (3.5) 

= 2{z + l)k' + {z + 1)/' + 2(z - 1)/ + {z- l){2z - 4)j. (3.6) 

The first two of these equations are second order, while the third one is first order. This 
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means that there must be five integration constants: Ci, ..., C5. The solution is given by 

j(r) = _ (i±l)^ g-(^+2)r _ (^+ ^)g2 ^_i(^+2+/3,)r ^ (^ + 1)^3 ^_ i (2+2-/3,)r^ j^g^^^j 

f rr^ - ^^^ p-(^+2)r , 2(5Z - 2 - f3,)c2 i(^+2+;3,)r _ 2(52: - 2 + f3,)c3 i(^+2-/3,)r , „ 

^^ ^ z + 2 ^ z + 2 + f3, z + 2-f3, + "' 

k(r) = ^^ e^(^+2)^ - 2(32:-4-/3,)c2 i(.+2+/3.)r , 2(3z - 4 + /3,)C3 i(,+2-/3.)r. 

^^z + 2 ^ + 2 + /3, 2 + 2-/3, 

This solution holds for the range of the dynamical exponent 1 < z < 2 and z > 2. The 
case of 2 = 2 must be treated separately because of the appearance of logarithmic modes. 
To see that these logarithmic modes are needed to solve the field equations, notice that 
e.g. the Ci and C2 modes have the same radial behavior when z = 2 (and similarly for the 
C3 and C4,5 modes). The case of 2; = 2 shall be discussed later on, in section 3.5. It is also 
interesting to see that the vector modes decouple from the metric modes when z = 1 (just 
rescale ci,2,3 -^ {z — l)ci,2,3 and shift C4 and C5). 

For 1 < 2; < 2, all the modes decay as r — ;■ 00; when z > 2 the C3 mode diverges. Notice in 
particular that C4 and C5 correspond to linearized diffeomorphisms, generated by the vector 
field 

e = |t9, + |x^9,. (3.8) 

Since a is a scalar, it cannot depend on C4 and C5 at linear order. In fact it is easy to see 
that, at this order 

a-ao = ej{r) + 0{e^). (3.9) 

The ci mode should be related to the mass of a black hole solution, as suggested by the 
asymptotically Lifshitz black holes considered for example in [14,21]. Additional evidence 
for this is provided by the ADM mass: 

4:(z — 2) 
Madm = e ^ ,' ci. (3.10) 

Notice however that this computation is a bit suspicious because the ADM mass is com- 
puted by means of background subtraction, and additional counterterms could modify the 
answer. 
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Second-order solution. The second order field equations are given by 

= z(z + 2) (/ - 2]f + ^^r + iff + if + i^jf + S^jf + ^^^±^ ff 

+ 8,(, + 2)/_l£(i±s/+8./-i<i±^*'-^f-'-ii-^'' (3.11) 

z — I z — \ z — \ z — \ 

= ^(^ + 2)(/-2jy + r + 4//' + 4/2 + 4^^7' + 8^jy 

2; — 1 z — \ z —\ 

04.7 ~ 9,7 

+ 8z(z + 2)j+Az f + 8z]' + ^^ ~k' + —— ~k", (3.12) 

z — 1 z — \ 



z(z 



- 2) (/ - 2jf + f + 4]'f' + if + 4zjf' + Szjj' 



2; — 1 z — 1 z — \ z — \ 

,~ 4z(z + l) ~, ~, 8z(z + l)r, 

+ 8^(z - 2) J + ^_^ ^ /' + 8^/ + ^_^ ^ A;'. (3.13) 

Just like the first-order equations, these consist of one first-order differential equation and 
two second order ones, so again there are five integration constants. Another thing one can 
read off from these equations is that the only modes that can appear in the second-order 
functions j, / and k are products of the modes we had already found at first order. The 
solution is thus given by 

J(r) = Ji e-(^+^)^ + 22 e-^(^+^+^)^ + j3 e-^(^-+2-^)'- + J4 

+ J8 e-5(3(^+2)+/3> ^ ^-^ ^- 1 (3{2+2)-/3)r^ 

/(r) = /i e-(^+2)r ^ ^^ g-i(.+2+/3)r ^ ^^ g-i(.+2-/3)r- ^ ^^ 

+ /5 e-(^+2+/3)r ^ ^^ g-(.+2-/3)r ^ ^^ g-2(.+2)r- ^g^^g) 

^(r) = fci e-(^+2)^ + ^2 e-^(^+2+'^)^ + /ts e-^(^+2"'^)'- + k^ 

+ ^5 e-(^+2+/3)r ^ ^^ g-(.+2-/3)r. ^ ^^ g-2(.+2)r ^3_ ^g) 

The coefficients jj, /«, and /cj depend on the dynamical exponent z as well as the first and 
second-order integration constants Cj and Cj. Instead of listing the coefficients explicitly, we 
shall discuss which coefficients are fully determined by the field equations and which ones 
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are related by integration constants. The coefficients of the modes that did not appear at 
ffist order are entirely determined by the field equations and thus only depend on the first- 
order integration constants q.^ The coefficients (j,, /,, ki) are related by q, with i = 1,2, 3. 
The coefficients f^ and /C4 do not enter the field equations at all, so let us call them f^ = C4 
and k4 = C5 for consistency of notation. 



3.2 On-shell action 

These solutions allow us to compute the on-shell action as a function of the integration 
constants up to second order in the expansion parameter e, 

r = ^(o)+£5(i)+£2^(2) + ... (3.17) 

The leading-order term vanishes, 5'(o) = 0, while the first-order term is given by 

which reproduces the result of [12]. At second order in e we find 

2ciC4(^-l) , 2cic,z{z + l) , 2c2C3(z + l)((z-2)/3-8) ^ , 

^(^) = 7T^ + — 7T~2 — + (.-2)(^-i) + "^' ^^-^^^ 

where ci is a second-order correction to ci, which cannot be determined by the asymptotic 
analysis and is therefore arbitrary. It is pleasing to see that the second order on-shell action 
is finite. Furthermore we recognize a familiar structure source/state, where C3 sources C2 
while a linear combination of C4 and C5 sources Ci. Therefore we identify C3, C4 and C5 
with the boundary conditions in the constant scalar perturbations sector. The remaining 
parameters ci and C2 are then naturally identified with the state of the system, and are 
determined by the initial and final conditions [22]. 

Since the counterterms in [12] had been devised to cancel first-order divergences only, one 
does not expect them to properly renormalize the on-shell action at higher orders. At 
second order, one finds indeed that the on-shell action is infinite when one uses those 
counterterms. 



Third-order on-shell action. The second order solutions should be sufficient to check 
finiteness of the on-shell action at third order (we expect only third-order corrections to 
Ci coming from the third order solutions). In fact an explicit computation shows that the 
third-order contribution S'(3) is finite with our counterterms. This provides an additional 
non-trivial check that our counterterms indeed remove the divergences. 



^The coefficients we are talking about here are J4 and (jj, fi, ki) with i = 5, ..., 9. 
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3.3 Vector and tensor perturbations 

For vector and tensor perturbations, the linearized analysis should be sufficient to compute 
the second-order on-shell action. Using again the notation in [12], we have: 

7t. = -e'^'^i.(r) + e%2i(r), (3.20) 

7., = e2^fc,,(r), (3.21) 

A, = y3^e^^yH(r), (3.22) 

where 



tdir) toir) 
to{r) -U{r) 



%(0 = ( Tli r,L 1 • (3-23) 



Vector perturbations. The vector sector is parametrized by 

vu{r) = cu + C2.e-(^+2)r ^ ^^^^-s.r-^ ^3 24) 

V2:{r) = 4^^C2.e(^-^)^ + ^c3.e-(^+2)^ + c,,. (3.25) 

z[z — 4) 2 + 2 

The on-shell action (3.17) converges only for z < 4, and its second-order (in e) term is 
given by^ 

4(z -1)(^ + 1) 2(z-l)(z^-4 z 

[7+2^ '''''''' + z{z-A) 

Once again we recognize the source/state structure cu/c^i and C4j/c2i. 



o(2) = ^_, , ^^ ciiCsj H ^j- Y\ C2iC4j. (3.26) 



Tensor perturbations. Finally the tensor modes, given by: 

td(r) = t,i + td2e-(^+2)^ (3.27) 

to(r)=tol+to2e^^'+'^^ (3.28) 

lead to the second-order contribution to the on-shell action (3.17) 

^(2) = 2(Z + l){tolto2 + tdltd2). (3.29) 

The source/state structure is in this case tdi/td2 and toi/to2- 



3.4 Boundary conditions 

At this point we can draw some conclusions: the modes C4, C5, Cij, C4J, t^i and toi should be 
interpreted as the sources in the metric sector, because they change the boundary values 



'For simplicity, we take the perturbation to lie along the x— axis. 
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of 7jj. The expectation values of the dual operators are given by Ci, c^i, C2i, td2 and to2- 
The mode C3 is a source for the massive vector field and is fixed by the leading term in 
a — oq. The corresponding expectation value is given by C2- 

We note that the mode C2i is problematic when z > 4, since it leads to a divergent on-shell 
action. 



3.5 The special case of 2; =2 

Let us repeat the analysis of the constant (scalar) perturbations for z = 2. The first-order 
solution was computed in [12]; it is given by 

j(r) = - (ci + C2 r) e"^'' + C3, 

/(r) = ^ (4ci - 5C2 + 4c2 r) e"^'' + {Ac^ r + C4) , (3.30) 

k{r) = — (4ci + 5c2 + 4c2 r) e"^'' + (-2c3 r + C5) . 

In this case, the modes ci and C2 are normalizable, while C3, C4, and C5 are non-normalizable. 
The possible modes in the second-order solution can again be obtained by squaring the 
first-order modes. The second-order solution is thus given by 

j{r) = [Ji+j2r + j3r ) e + Ji+ J5r + Jer + [jj + jsr + j^r ) e , 

/>) = (/i + /2r + hr') e-'^ + f, + f^r + f^r' + {fj + fsv + hr^) e-«^ (3.31) 

k{r) = {ki + k2r + k^, r^) e"^'' + ^4 + /cs r + /cg r^ + (/cy + fcg r + /cg r^) e"®''. 

Again, the coefficients ji, fi, ki depend on the first-order and second-order integration 
constants Cj and Cj respectively. We shall not list these coefficients explicitly, but let 
us mention where the second-order integration constants appear. The coefficients of the 
modes that were not present in the first-order solution are all fully determined by the field 
equations.^ The coefficients {ji,fi,ki) are related by the integration constant Ci, while 
iJ2, /2, k2) are related by C2, and (j4, f^, k^) are related by C3. The coefficients f\ and k^ do 
not enter the field equations at all, so just like the z ^ 2 case, we call them /4 = £4 and 
A;4 = C5 for consistency of notation. 



The on-shell action (3.17) for z = 2 at first order in e is given by 

5(1) = ^, (3.32) 

which reproduces the result in [12]. At second order, we find 

25 1 25 19 , , 

S(2) = 3CiC3 + — C2C3 - -C1C4 - 777C2C4 + C1C5 + —:C2C5 + C2, (3.33) 

12 D 24 12 



To clarify, these coefficients are j^ and {ji, fi,ki) with i = 3, 6, 7, 8, 9. 
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where C2 is a correction to the vev C2. We also checked whether the on-shell action is finite 
at third order and we find that it is. 

Boundary conditions The boundary conditions are essentially the same as in the z ^ 2 
case. The modes C3, C4, C5 should be seen as the sources, while the vevs are represented by 
Ci and C2. Again, C3 can be interpreted as the source for the vector's mass term via (3.9). 
Notice that the C2 mode takes on a similar role as ci for -2 7^ 2, for instance the ADM mass 
is Madm = 4C2/3 when z = 2. 



4 Conclusions 



We have found a new and systematic method for simultaneously determining the bound- 
ary conditions on the one hand, and finding the counterterm action for asymptotically 
Lifshitz spacetimes on the other hand. This method allowed us to find contributions to 
the 'Lifshitz scaling anomaly'. We performed a non-trivial consistency check for the coun- 
terterms obtained via this method. We saw that our counterterms properly renormalize 
the on-shell action even for higher-order perturbations. The counterterm action we find is 
a local functional of the fields by construction. Moreover, we find that the counterterms 
are independent of the radial cut-off, unlike some previous approaches, e.g. [23] for z = 2. 
Although in that paper the counterterms are local functionals of the boundary fields, the 
coefficients that appear depend explicitly on the radial cut-off, whereas in our case this 
cut-off dependence is only implicit through the dependence of the fields on the radius. The 
counterterms found in [12] only managed to make the on-shell action finite up to linear 
order. In our case, the on-shell action remains finite when we turn on non-normalizable 
modes. 

It should perhaps be emphasized that the idea that one should be able to remove all 
divergences using only local counterterms is a conjecture and we have only shown that it 
works in particular examples. A full and general proof that this works for all reasonable 
bulk Lagrangians and to all orders is lacking. 

Although we focused on constant perturbations, the Hamilton- Jacobi analysis can be used 
to find higher-derivative counterterms as well. See Appendix B for an example of such a 
calculation. Though our method determines the asymptotic behavior of the fields, there 
are still some puzzles that remain. For example, in the computations in section 3, five free 
parameters appeared, whereas for a non-degenerate set of equations of motion one would 
expect to find an even number that can be split in "coordinates plus momenta" or equiv- 
alently in sources and expectation values, as was emphasized in [24]. One would expect 
that a canonical analysis in this sector would reveal that one of the five parameters can be 
removed by a suitable gauge transformation which seems related to a bulk diffeomorphism 
(which turns out to be a Lifshitz rescaling at the linear level). However, a preliminary anal- 
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ysis suggests that the corresponding constraints will end up being non-linear once higher 
order corrections are included, and the precise nature of these non-linear boundary condi- 
tions remains to be determined. It is also unclear whether this constraint would somehow 
follow from the analysis of the Hamilton- Jacobi equations or require separate input. 

There are various qualitative differences for different values of the dynamical exponent z. 
As mentioned in section 3, for z > 4, divergences seem to appear which cannot be canceled 
using local counterterms. If we blindly follow the strategy we have been employing, we 
would be forced to impose more stringent boundary conditions for 2; > 4 which remove 
these divergences and it would be interesting to explore this in more detail. Furthermore, 
the source for the massive vector field is irrelevant when 2; > 2, and therefore we expect 
that non-local counterterms are needed at sufficiently high order in the sources to make 
the on-shell action finite, as pointed out in [25]. 

Another thing which would be interesting to compute is the holographic Weyl anomaly 
in a curved background for the Lifshitz case, which from a preliminary analysis seems to 
involve terms like RabA°'A^ for z = 2 and d = 3. Though in principle straightforward, the 
relevant computations turn out to be extremely tedious and we leave this for future work. 

One may wonder whether the fact that all divergences can be canceled by local counterterms 
is a special feature of field theories with a Lifshitz dual or hold for a more general class of 
non-relativistic scale-invariant theories, and it would be interesting to explore this question 
directly in field theory. 

There are many further directions to explore, such as applications to black hole solutions, 
applications to correlation functions, the extension of our work to Schrodinger spacetimes, 
etc, and we hope to turn back to some of these in due course. 
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A The Einstein— Proca Hamiltonian 



In this section we shall compute the Hamiltonian associated to the Einstein-Proca action, 
which is given in section 2.1. We denote by S,. the surface of constant radial coordinate 
r. We assume for simplicity that dUr = 0, so that we need not worry about possible 
boundary terms later on. The foliation can be written in the form of a parametric relation 
X^ = X^{r,x°'). It is useful to define the projector 

dXi" 

which projects onto the directions tangent to the hypersurface. Thus, the projector is 
orthogonal to the unit normal, n^p^ = 0. The cotangent basis is spanned by 

dX'' = r''dr + p^dx''. (A.2) 

The vector r'* points along the radial flow, which does not necessarily mean it should be 
proportional to the unit normal. The flow vector is generically given by 

r^" = Nnf" + N"" p^, (A.3) 

where the normal and tangent pieces are given by the lapse A^ and shift A^" respectively. 
The metric Q/j,!, can be rewritten in terms of the fields (A^, A^", '-fab) as follows. 

ds^ = N^ dr^ + -fabiN" dr + dx''){N^ dr + dx^), (A.4) 

where 7afe = PaPb di^i' is the induced metric or first fundamental form. In order to rewrite 
the gravitational Lagrangian in terms of quantities that are either intrinsic or extrinsic to 
Sr, we use a projected form of the Gauss-Codazzi equations, namely 

^(d+i) ^ R + K^-KabK''^ + 2V^{rfVun^'-n^'Vun''). (A.5) 

The extrinsic curvature, or second fundamental form, is given by Kab = \.^nlab = VaVl V/^ni, 
Similarly, the Maxwell term and the mass term can be split up into normal and tangential 
pieces using the completeness relation g^'^ = n^n'^ + PaPblab- Let us use the short-hand 
notation V = n'^A^ and K-a = -^nAa — daV, such that 

F^^F^" = FabF'''' + 2KalC\ (A.6) 

A^A^" = AaA'^ + V^. (A.7) 

The actions from before thus become^ 

^-grav = 1^ [dr f d^Xy/^N [R-2k + K^- KabK'"") , (A.8) 

Sa = j dr j d'^x,/^ N(--^FabF'^' ~]^m^ AaA- -]^lCalC'' -]^m^vA. (A.9) 



^By the Lie derivative of a tangential object with respect to some (d + l)-dimensional vector ^'', e.g. 
■^?^ai...a„, we reaUy mean p^J ■■■PaZ •^«^pi-a'>.- 

^We take E^ to be the boundary of our {d + l)-dimensional space such that the total-divergence term 
in (A.5) precisely cancels against the Gibbons-Hawking term in the gravitational action. 
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Let us define the Lagrangians such that S'grav = J dr Lg^^^ and Sa = J dr La, via which 
we may obtain the Hamiltonians by means of a Legendre transformation. Before we do so, 
however, we must define our generahzed velocities first. 



Aa = ^rAa = NJCa + N'F.a + daiNV + N'At) 



Finally, the canonical momenta are^ 

1 SL 



ab 



IT 



E" 



grav 



1 SLa 



2k^ 



{K^'-K-f 



ab\ 



-JC". 



V -7 5Aa 
Now, we are ready to perform the Legendre transformation 



(A.IO) 
(A.ll) 



(A.12) 
(A. 13) 



H, 



;rav - / d'^ X y/^ Tl"^ ^ ab - L ^ 



d'^x^^ { N 



0,-2 I ^ ^ab 



d- 1 



TT 



^(^-2A) 



2k 



(A.14) 

Ar"(-2DV„,) 
(A.15) 



and similarly for the vector field 



H, 



d'^x^/^E'^Aa-L, 



(A.16) 



d'^x^^ { N 



ra 



^EaE'^ + ^FatF'^' + Im^ AaA'^ + ^ ^' " ^D^E'^ 



+ N'^{FatE'-AaD,E') 

(A.17) 



We can also combine the above two Hamiltonians as 



H = Hg,^, + HA = I d^xy^^^NH + N-Ha^ 



(A.18) 



where we introduced the Hamiltonian constraint "H and the momentum constraint Ha- 
After integrating out the non-dynamical field V, the Hamiltonian constraint function is 



n 



0,-2 / ^ ^ab 



d- 1 



-TT 



^EaE'^ - ^{DaE'^f - C. (A.19) 



^Note that, strictly speaking, the canonical momenta are -y/— T'tt"'' and y/—jE°-. This will be taken 
into account in the Legendre transformations to be performed in (A.14) and (A.16). 
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Here, L is the Lagrangian density restricted to the hypersurface S 



ri 



C = -^(i? - 2A) - -KbF'^' - \m' A^A\ (A.20) 

The momentum constraint function is given by 

Ua = -2D^'Ka, + Fa,E'-AaD,E\ (A.21) 



B Higher-derivative Counterterms 

In this section, we briefly mention how one could systematically solve the local part of 
the Hamiltonian constraint at the level of higher derivatives (thus finding the local higher- 
derivative counterterms). We shall put the scalar field = 0, as it will not be more 
illuminating in this specific discussion. 

Local on-shell action Ansatz. In the following we are interested in deformations that 
involve only the metric 7ab and the massive vector Aa. 

Aoc = U{a) + C{a)DaA'' + P(«)AM''D,A + ^{a)R + . . . (B.l) 

Of course, there are other two-derivative terms as well as higher-derivative terms in the 
Ansatz, but for our purpose of illustrating our method these terms will suffice. We assume 
that (9Sr = 0, so that we need to specify the possible counterterms only up to total 
derivatives. We perform a derivative expansion, 

4o°c = U{a\ (B.2) 

£2 = C{a)DaA'' + V{a)A''A'DaAk, (B.3) 

4oc = *(«)^+--- (B-4) 

and 

£(0) = _2A _ !^«, (B.5) 

C^'''^=R-\FabF-\ (B.6) 

The non- derivative level (level zero) has already been covered in Section 2, so let us go 
directly to the level of one spacetime derivative. 

One derivative. At level one we have only two possible structures. The canonical mo- 
menta are given by 

.(1) 



1 'i^ 



-V - C J {A''A\D ■ A) - 2A'A^''D''^A^ + -f''\A''A'^D^Ad)) 



-7 7ab 

(B.7) 
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and 

^(i)a _ (2C' - V) {A^iD ■ A) - AbD'^A'') . (B.8) 

The Hamilton constraint can be solved if 

V = 2C'. (B.9) 

The resulting term in C\J^ is just a total derivative, 

CDaA"" + 2C'A''A''DaAb = DaiCA''), (B.IO) 

and can be discarded. 

Two derivatives: the $i?, term. Since the $/? term does not mix with the other 
two-derivative terms, we can consistently solve for 0(a). 

-'='■" = >42-f^ + ... (B.n) 

^(2)a ^ 2^'RA'' + ... (B.12) 

We now want to compute the coefficient of the R term. Only terms with i? or a not 
contracted R^y can produce a R term in the final expression: 

^1? = \9ab^R - ^'AaAbR - Rab^ + . . . , (B.13) 

^(2)a ^ 2$'i?A" + . . . (B.14) 

Therefore we have 

2{S[^^,S[2} - C^^^ = R i-j^U + -A^{^'U - ^U') + (4^2 + A^)^'U' + l") + . . . = 0. 

(B.15) 
Again, we expand $ in power series in (a — oq) where a = A?, 

$ = 6o + 6i(a - tto) + &2(« - «o)^ + • • • , (B.16) 

and we plug this result into (B.15). We obtain 

ho = \. (B.17) 

A similar computation for hi yields 

h - 5Z-2 + /3, 

^^ - 4(. + l)(.-2 + /3.)- ^^-^^^ 

There does not seem to be a continuous ambiguity for the higher order coefficients. 
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Let us briefly discuss an important feature of (B.15). The function $ satisfies a first-order 
differential equation, tlierefore it seems somewhat strange that we were able to determine 
the coefficient bo uniquely, which amounts to specifying the initial condition. The reason 
for this is that, since we want to compute the polynomial part of the on-shell action, we 
are using a power-series expansion. Nevertheless, the general solution of the differential 
equation might not be polynomial, so by requiring that our solution is a polynomial, we 
are effectively determining the initial condition. We illustrate this phenomenon with a toy 
example. Consider the differential equation: 

xf'{x) + af{x) + 1 = 0. (B.19) 

If a 7^ this has the following general solution: 

-- + Ax"", (B.20) 

(Jj 

where A is an arbitrary constant. If a 7^ 0, —1, —2, . . ., then the solution is not polynomial 
and using a Taylor expansion amounts to choosing A = 0. Nevertheless, if the coefficient 
a is a negative integer, the solution is indeed a polynomial but A is undetermined. This 
amounts precisely to a continuous ambiguity that we would find by using the power-series 
method. 

Equation (B.15) can be cast in a form similar to the toy model we just considered: 

a^ + 4a)f/' + ^aUj $' + (^aW - jUj $ + 1 = 0. (B.21) 

The coefficient of $' is simply —dra, and stability required that this coefficient vanishes 
as a — )■ ttO) ^^ we explained at the end of section 2.3. This feature is very general and it 
explains why the HJ method is able to fix the derivative counterterms. 
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